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Abstract 


The cloud radiative effect (CRE) of each longwave (LW) absorption band of a GCM's 
radiation code is uniquely valuable for GCM evaluation because (1) comparing band-by-band 
CRE avoids the compensating biases in the broadband CRE comparison and (2) the fractional 
contribution of each band to the LW broadband CRE (f C RE) is sensitive to cloud top height but 
largely insensitive to cloud fraction, presenting thus a diagnostic metric to separate the two 
macroscopic properties of clouds. Recent studies led by the first author have established 
methods to derive such band-by-band quantities from collocated AIRS and CERES observations. 
We present here a study that compares the observed band-by-band CRE over the tropical 
oceans with those simulated by three different atmospheric GCMs (GFDL AM2, NASA GEOS-5, 
and CCCma CanAM4) forced by observed SST. The models agree with observation on the 
annual-mean LW broadband CRE over the tropical oceans within ±lWm" 2 . However, the 
differences among these three GCMs in some bands can be as large as or even larger than 
±lWm" 2 . Observed seasonal cycles of f C RE in major bands are shown to be consistent with the 
seasonal cycle of cloud top pressure for both the amplitude and the phase. However, while the 
three simulated seasonal cycles of f C RE agree with observations on the phase, the amplitudes 
are underestimated. Simulated interannual anomalies from GFDL AM2 and CCCma CanAM4 are 
in phase with observed anomalies. The spatial distribution of f C RE highlights the discrepancies 
between models and observation over the low-cloud regions and the compensating biases from 
different bands. 
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1. Introduction 


Since the 1980s, broadband radiative flux and cloud radiative effect (CRE, the difference 
between all-sky and clear-sky fluxes), have been extensively used in climate studies (e.g. 
Ramanathan et al. 1989; Wielicki et al. 2002; Wong et al. 2006), especially in the evaluation of 
climate models and cloud feedback studies (e.g. Allan et al. 2004; Allan & Ringer 2003; Raval et 
al. 1994; Slingo et al. 1998; Wielicki et al. 2002; Yang et al. 1999). In the development of a GCM 
for climate studies, one inevitable and important step is "tuning" in which poorly constrained 
parameters are adjusted using observations and physical principles to ensure energy balance at 
the top of atmosphere (TOA). This ensures that simulated broadband quantities are generally 
consistent with the observed counterparts at the TOA. However, it cannot guarantee the 
consistency between band-by-band decompositions of observed and simulated radiation fluxes. 
In fact quite often the compensating biases from different absorption bands offset each other 
and lead to a seemingly good agreement between modeled and observed broadband fluxes 
(e.g., see the work in the thermal infrared by Huang et al. 2006; Huang et al. 2007; Huang et al. 
2008; and Huang et al. 2010). Similar compensating biases can be expected in the simulated 
broadband CRE as well. Therefore, directly using of band-by-band flux and CRE in model 
evaluation can avoid the compensating errors and highlight the biases in different bands since 
flux and CRE of each individual molecular absorption band are calculated directly in the GCM 
radiation scheme. Moreover, as illustrated in Huang et al. (2010) with both conceptual model 
and numerical simulation, the fractional contribution of each band to the broadband longwave 
CRE (for brevity, hereafter denoted as f CRE ) is sensitive to cloud top height but largely insensitive 
to cloud fraction. The LW broadband CRE can be written as CRE L w = fc (F c i r - F c id) where F is the 
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outgoing longwave flux at top of atmosphere, fc is the cloud fraction, and subscripts c ir and c id 


denote clear-sky and overcast cloudy sky, respecitvely. Correspondingly, the i-th band CRE can 
be written as CRE 1 = fc (F' C | r - F' C | d ), and the superscript 1 denotes the i-th band. The resulting 
expression for/ CR f is. 


fc 


CRE 


F' dr —F‘cld 
^clr ~ F C ld 


( 1 ) 


Therefore, the common factor of cloud fraction cancels and f C RE is only sensitive to cloud top 
temperature, making f C RE a useful quantity in diagnosing and evaluating modeled CRE. The LW 
broadband CRE is sensitive to both the cloud fraction and (mostly) cloud top height while the 
SW broadband CRE is sensitive to both the cloud fraction and (mostly) the cloud water path (i.e. 
cloud reflectance). Since the LW f C RE is sensitive to cloud top height but not to cloud fraction, it 
provides a dimension to sort out the contributions of cloud faction and cloud top height to the 
broadband CRE. 


Huang et al. (2008) and Huang et al. (2010) mainly focused on the algorithm for deriving 
such band-by-band LW flux and CRE from the AIRS spectral radiances collocated with CERES 1 
observations from the same Aqua spacecraft. Using one year of data and corresponding 
simulations from a GFDL AGCM (AM2), both studies also showed preliminary usage of such data 
in GCM evaluation. To further explore and demonstrate the potential of such band-by-band CRE 
in model evaluations, this study performs the first ever comprehensive evaluations of GCM- 
simulated band-by-band CRE from three different GCMs. Multiple years of data are used to 


1 AIRS stands for Atmospheric Infrared Sounder, which is a grating spectrometer. CERES stands for Clouds and 
Earth's Radiant Energy System, which consists of two broadband radiometers and one narrow-band radiometer. 
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explore the long-term mean, seasonal cycles, as well as interannual variations of CRE and f C RE of 


each LW band. 

The data and GCMs used in this study are described in Section 2. Sections 3 and 4 
present the comparison results. Conclusion and further discussion are given in Section 5. 

2. Observations and GCM simulations 

2.1 Observations 

The algorithms in Huang et al. (2008) and Huang et al. (2010) were developed for and 
validated against observations over the tropical open oceans (30°S to 30°N). Hence here we 
employ the algorithms to derive clear-sky and all-sky spectral fluxes at a 10cm" 1 interval for the 
longwave spectrum (0-2000 cm" 1 ) using the collocated AIRS and Aqua-CERES observations over 
the tropical open oceans from 2003 to 2007. The monthly mean spectral CRE is then calculated 
in a similar way as the broadband CRE derived from ERBE or CERES observations. The scene 
type information from CERES SSF (Single Satellite Footprint) is used with pre-developed spectral 
anisotropic distribution model (ADM) to invert spectral flux for each AIRS channel, and a multi- 
linear regression scheme is then used to estimate the spectral flux over spectral region not 
covered by the AIRS channels. As shown in Huang et al. (2008) and Huang et al. (2010), the OLR 
(outgoing longwave radiation) derived by this approach agrees well with the collocated CERES 
OLR and this good agreement is found over different cloud scene types (as distinguished by 
distinct cloud fractions and cloud-surface temperature contrasts). Comparisons with synthetic 
data showed that the algorithm can reliably estimate spectral flux at 10cm" 1 resolution with 
maximum fractional difference less than ±5% for clear-sky scenes and ~±3.6% for cloudy-sky 
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scenes. Details of the algorithm and validation can be found in Huang et al. (2008) and Huang et 


al. (2010). 

Table 1 summarizes the comparison between OLR derived by this spectral ADM 
approach with the collocated CERES OLR for 2003 to 2007. AIRS version5 calibrated radiances 
and CERES SSF (single satellite footprint) Edition 2 data are used. The 2a radiometric calibration 
uncertainty of CERES OLR is ~1%, which translates to ~ 2.5Wm' 2 for a typical OLR value in the 
tropics. Tablel shows that, for all the years, the mean difference between CERES clear-sky OLR 
and our derived clear-sky OLR is always below the 2a uncertainty. The cloudy-sky OLR mean 
difference is about the same as the 2a uncertainty. The standard deviations of these 
differences change little from year to year. 

2.2 GCMs 

We used three atmospheric general circulation models (hereafter, GCM) in this study, 
the GFDL AM2 model (GFDL GAMDT 2004), the Fortuna 2.2 version of the NASA GEOS-5 model 
(Molod et al. 2012), and the fourth generation AGCM at CCCma, Environment Canada (CanAM4) 
(vonSalzen et. al. 2012, manuscript submitted to Atmosphere-Ocean). Each GCM is forced with 
observed SSTs over multiple years. For each model, clear-sky flux and all-sky fluxes of the 
individual bands of the LW radiative transfer code are directly saved in the output of 
simulations. For GEOS-5 and CanAM4, the simulations were carried out for 2000-2009 while 
only the 2003-2007 period is analyzed and compared with observations here. For AM2, due to 
the limited resources, the simulation was carried out only up to 2005 so only the simulations 
from 2003 to 2005 are analyzed here. As will be shown in the next section, the year-to-year 
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variations of band-by-band CRE and f C RE are small so the different period in the AM2 


simulations has little effect on the comparisons of long-term mean and mean seasonal cycle. 

The longwave radiation scheme in the GFDL AM2 model follows Schwarzkopf & 
Ramaswamy (1999). The LW spectrum is divided into eight bands, with two water bands (far-IR 
and >1400cm _1 ) treated together, in practice. Clouds are assumed to be non-scattering in the 
LW. The LW radiation scheme in the NASA GEOS-5 is based on Chou et al. (2003). It divides the 
LW spectrum into 9 bands and can be run at two accuracy modes (high accuracy mode is used 
here) using either a k-distribution or look-up table approach to calculate gaseous transmission 
functions. Cloud scattering in the LW is handled empirically using a rescaling approach. In the 
CanAM4 model, the optical properties of gases is modeled using a correlated-k distribution (Li 
and Barker 2005) while the radiative transfer is simulated using the Monte Carlo Independent 
Column Approximation and two-stream radiative transfer solutions (von Salzen et al. 2012, 
manuscript submitted to Atmosphere-Ocean; Pincus et. al. 2003; Li 2002). Scattering by clouds 
droplets is included based on the absorption approximation (Li and Fu, 2000). Note both GEOS- 
5 and CanAM4 models take into account the scattering of cloud in the longwave. 

2.3 Redefining the bands for comparison 

Because of the different radiation schemes used in the three GCMs, the number of 
bands and the bandwidth for each band are not necessarily the same. In fact, GFDL AM2 uses 8 
bands in the LW while NASA GEOS-5 used 10 bands and CanAM4 uses 9 (Table 2). To facilitate 
the comparison, we define five new bands (Table 2) to be used for comparisons between GCMs 
and between GCMs and observations, each of which is either a band common to the GCM 
radiative transfer schemes or some combination of bands used in the GCMS. This ensures the 
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maximum compatibility among GCMs for the band-by-band comparisons. As shown in Table 2, 


after merging into five bands, the bandwidth structures of GEOS-5 and CanAM4 models are for 
the most part consistent with each other and consistent with the five bands defined for the 
observational analysis. AM2 bandwidths are slightly different. To understand the effect of such 
intrinsic differences in the bandwidths in the GCMs on the band-by-band fluxes and CREs, we 
compute the flux of each band (F m od e i_bnd) for a blackbody at a temperature T using the 
bandwidth from each GCM. This flux is then compared to the flux of the band defined for this 
study (Fbnd)- The relative difference, (F m0 deij>nd-Fbnd)/Fbnd, is shown in Fig. 1 for T =220K and T= 
298K, respectively. Except for Band4 (ozone band, 980-llOOcm" 1 ) and Band5 (1100- 1400cm' 1 ), 
the relative difference due to the varying bandwidths is within 10% (Figurela and lb). By 
subtracting the flux for T=298K from the flux for T=220K, which mimics the magnitude of LW 
CRE (Fig. lc), we find that the relative difference for each band in each GCM's radiation scheme 
is about the same. The largest relative difference is ~ 33% and always occurs in the ozone band. 
This larger error occurs because of the rapid decrease of the blackbody curve between 220 and 
298K over this band, which causes the small bandwidth differences between AM2 and 
GEOS5/CanAM4 to lead to large fractional difference in flux. For this reason, we omit here the 
ozone band (band4) in some analyses presented in Sections 3 and 4. If an ozone-band result is 
discussed, it must be interpreted with caution for the GFDL AM2 which has a different 
bandwidth for the band from the other two GCMs and the observations. For most of the 
analysis, our focus is primarily on bandl, band2, and band3, for which the intrinsic bandwidth 
differences in the GCMs yield no more than 10% relative difference in the flux, CRE, and f CRE . 
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3. Model intercomparisons and comparisons with observations: averages over the tropical 


ocean 

Figure 2 shows the five-year averages of all-sky OLR and LW CRE from the NASA GEOS-5 
and CCCma CanAM4 model simulations, three-year averages from the GFDL AM2 simulation, 
and the counterpart 5-year averages from the collocated AIRS and CERES observations 
(hereafter, referred to as observations for the brevity). Generally the models and the 
observations are consistent with each other on the broad spatial features of both all-sky OLR 
and longwave CRE. In the tropical Pacific warm pool east of the Maritime Continent, the AM2 
and CanAM4 models have lower all-sky OLR than the observation and the GEOS-5 model. 
Accordingly, their LW CREs are higher than those in the observations and the GEOS-5 model. 
The LW CRE is only ~10 Wm" 2 (or even smaller) over the oceans west of major continents 
because these regions are frequently covered by marine stratus with cloud top ~lkm and thus 
has little thermal contrast with surface. Judging from the LW broadband CRE here alone, 
observations and all GCMs seemingly agrees better over the low-cloud regions than over the 
high-cloud regions such as the ITCZ and SPCZ (Southern Pacific Convergence Zone). However, as 
we shall show in the following subsections, the f C RE over the low-cloud regions is indeed quite 
different among the GCMs and between the models and observation. 

3.1 Long-term means and seasonal cycles 

Table 3 summarizes the multi-year means of the LW broadband and band-by-band CREs 
averaged over the entire tropical oceans. Observed long-term mean LW broadband CRE is 28.5 
Wm" 2 while the model results are slightly lower by 0.2-1.2Wm" 2 . However, the difference in the 
CRE of a particular band could be easily as large as, or even larger than, the broadband 
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difference. For example, the differences between observed and simulated CREs in Band2 (the 
C0 2 band) are ~-0.6 to 0.7 Wm" 2 . The difference between observed and CanAM4-simulated CRE 
in Bandl alone is ~1.12 Wm" 2 . Modeled broadband CREs are always smaller than the observed 
CRE, but modeled CREs in a given band could be either larger or smaller than its observed 
counterpart. Such compensation among bands leads to an apparent good agreement of the 
broadband LW CRE. Band-by-band CREs exposes such compensation in a quantitative way and 
makes it possible to further examine the sources of the compensating biases. 

Among the five bands examined here, the largest contributor to LW CRE is Band3, the 
first window region from 800 to 980 cm" 1 , and is responsible for ~32-36% of LW CRE. There is a 
large discrepancy in the CRE of Band4 (ozone band) between the AM2 model and the other two 
models and the observations. As mentioned in Section 2, the AM2 model has a different 
bandwidth for Band4 while others have identical bandwidth. Second, the AM2 simulation here 
is done with prescribed climatological ozone profiles at the 1990s level obtained from a 
combination of ozonesonde and satellite measurements (Paul et al. 1998; GFDL GAMDT 2004). 
In contrast, the GEOS-5 model ozone fields are calculated on-line inside the GCM using a 
parameterization described in Rienecker et al. (2008). The CanAM4 simulation used a zonally 
averaged version of AC&C/SPARC ozone database that was prepared for the CMIP5 (Coupled- 
model inter-comparison project, phase 5) model simulations (Cionni et al. 2011). Over the 
historical period, the database consists of time-varying ozone fields based on observations for 
the stratosphere and chemistry-climate-model simulations for the troposphere. Band 4, the 
ozone band, never saturates, which makes it sensitive to stratospheric and tropospheric ozone 
as well as to surface and cloud temperatures. Moreover, ozone spatial distribution is indeed 
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affected by the large-scale circulation and transport. Therefore, besides the band-width 


discrepancies, whether realistic on-line ozone fields are available or not could potentially affect 
the simulated CRE results as well. 

Seasonal cycles of CREs of Bands 1-3 are shown in Fig. 3. Compared to the mean values 
shown in Table 3, the seasonal fluctuation is very small (~3.5% or even smaller). For both 
models and observations, each band CRE seasonal cycle closely tracks one another. The 
observations show two peaks with one in April-May and the other in November-December, 
which is related to the movement of the Sun as well as the north-south seasonal movement of 
ITCZ. The models also show such semi-seasonal cycles with similar phase except that the first 
peak in the GEOS-5 is one month ahead of the observations. Although the phases are generally 
consistent to each other, the magnitudes of the seasonal cycles are noticeably different by 
about a factor of two, with the observations and GEOS-5 having similar amplitudes, while those 
of AM2 and CanAM4 being much smaller. 

The good phase syncing between CREs of different bands is largely due to the fact that 
for all the LW bands the surface and cloud top temperature contrast is the largest driver of CRE. 
Therefore, they all vary in phase with the LW broadband CRE in terms of the absolute 
magnitude. Seasonal cycle of the f C RE (Fig. 4) is different because, in this way, the LW broadband 
CRE is always normalized to 100% for each month. In contrast to the seasonal cycle of absolute 
band CRE (Fig. 3), for both models and observations a single annual cycle is dominant instead of 
a semi-annual cycle. The seasonal cycle of Bandl (water-vapor band) f CRE is strongly anti- 
correlated with that of Band3 (the window region) f C RE- For both models and the observation, 
Bandl peaks in March-April while Band3 peaks in July-August. This suggests that any phase 
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discrepancies shown in Fig. 3 are likely dominated by differences in cloud fraction rather than 


differences in cloud top height (otherwise, the same discrepancies would have appeared in 
such f CRE seasonal cycles as well). For f CRE seasonal cycles, the GCMs have similar amplitudes but 
are 2-3 times smaller than the observed amplitude (Fig. 4). 

Figure 4 also shows that the seasonal cycle of Band2 f CRE closely tracks that of Bandl but 
anti-correlated with that of Band3. Figure 5 plots the correlation coefficients between the 
seasonal cycle of Bandl f CRE and the seasonal cycle of f CRE of the remaining bands. Models and 
observations behave similarly in terms of such inter-band correlations: window bands (Band3 
and Band5) are strongly anti-correlated with Bandl. The observations exhibit a moderate 
positive correlation between Band2 and Bandl while much stronger positive correlations are 
seen for the GCMs. The largest discrepancies exist in Band4 (the ozone band). As expected, 
Band4 f CRE in observation, GEOS-5, and CanAM4 shows strong in-phase relation with Band3 f CRE , 
since cloud and surface thermal contrast significantly affects both bands. For the AM2 model, 
the correlation between Band4 and Band3 f CRE is weak, partly due to the prescribed 
climatological ozone profiles. Such differences indicate the benefit of having realistic or self- 
consistent ozone profiles in the transient simulations. 

3.2 Explanations of the seasonal cycles 

To understand the seasonal cycles shown in Figs. 3-4, it is instructive to look the 
seasonality of cloud macroscopic properties averaged over the tropical oceans. Fig. 6(a) shows 
the ISCCP (International Satellite Cloud Climatology Project, Rossow & Schiffer 1991; Rossow & 
Schiffer 1999) long-term mean seasonal cycles of cloud fractions for high, middle, and low cloud, 
respectively. The peaks of high-cloud fraction can be found during two periods: April-May and 
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November-January, the latter being also the peak for middle-cloud fraction. The seasonality of 


high-cloud and middle-cloud fraction is consistent with the observed seasonal cycles of band- 
by-band CRE in Fig. 3. Meanwhile, even though the low-cloud fraction peaks in July-August, it 
contributes little to the observed seasonal cycle of CRE because low clouds contribute little to 
the absolute LW CRE or the absolute band CRE. 

Fig. 6b shows the ISCCP mean seasonal cycle of cloud top pressure (CTP) averaged over 
the tropical oceans. In contrast to the seasonality of cloud fractions, it shows that mean cloud 
top is highest in April-May and lowest in August-September. As shown in Fig. 7b in Fluang et al. 
(2010), the higher the cloud top, the larger the Bandl contribution to the total LW CRE and the 
smaller the Band3 contribution. This behavior is mainly due to the shift of blackbody peak 
emission towards lower frequency as the cloud top temperature becomes colder (i.e. Wien's 
displacement law). Therefore, the seasonality of CTP is consistent with the observed seasonality 
of f CRE in Fig. 4. 

Given the strong anti-correlation between Bandl f C REand Band3 f C RE, it is worthwhile to 
explore to what extent we can use the ISCCP mean seasonal cycle of CTP to estimate the 
amplitudes of the seasonal cycles of Bandl and Band3. For both seasonal cycles, we use the 
standard deviation as a measure of the amplitude. Similar to Huang et al. (2010), we assume a 
layer of optically thick cloud (x»l) and typical tropical sounding profiles of temperature, 
humidity, and ozone (McClatchey 1972). Then we vary the cloud top from the lower 
troposphere to the upper troposphere and use MODTRAN5 (Berk et al., 2005) to calculate f CRE 
of Bandl and Band3 accordingly (Fig. 7a). When the cloud top is lower than 600 hPa, Bandl f CRE 
is less than ~6% while Band3 f CRE is more than 40%. As the cloud top moves upward, the 
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contribution from Bandl gradually increases while that from Band3 decreases. When the cloud 


top reaches 150 hPa, the contributions from both bands are essentially equal (fcRE~28%). Fig. 7a 
also shows that, above 600hPa, the change of the fractional contribution with respect to cloud 
top pressure is nearly linear. Therefore, we define the linear regression slope for the i-th band 
as 


_ dfcRE 
~ dCTP 


( 2 ) 


where superscript /' denotes the i-th band, ri and r 3 are then derived by linear regression 
of the curve shown in Fig. 7a. Then the amplitude of seasonal cycle of Bandl (Band3) can then 
be estimated by multiplying rl (r3) with the standard deviation of ISCCP CTP seasonal cycle in 
Fig. 6b, which is 12.2 hPa. Fig. 7b shows the amplitudes estimated in this manner vs. the real 
amplitudes derived from Fig. 4 top panel. The simple estimations are 82% and 85% of the 
observed values for Bandl and Band3, respectively. No cloud fraction values are involved in this 
estimation, which corroborates the fact that the fractional contribution of each band is largely 
insensitive to cloud fraction but is mainly affected by changes in cloud top pressure. This 
behavior is very different from that of the absolute amount of band-by-band CRE, which is 
largely influenced by the fractions of high cloud and middle clouds. 

Since the ISCCP CTP seasonality can largely explain the observed seasonality of f C RE, the 
model-observation discrepancies shown in Fig.4 should be largely due to the difference in 
observed and simulated seasonal cycles of CTP, especially the amplitudes. Using the same 
slopes of ri and r 3 derived above, we estimated the amplitudes (the standard deviations) of CTP 
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seasonal cycles for AM2, GEOS-5, and CanAM4 models are 2.9 hPa, 3.7 hPa, 4.7 hPa, 


respectively (in comparison to the 12.2 hPa derived from ISCCP's mean CTP seasonal cycle). 

3.3. Interannual variations 

Figure 8 shows the 5-month running means of deseasonalized time series of CRE along 
with the ENSO index. Following several previous studies (Deser & Wallace 1990; Fu et al. 1996; 
Klein et al. 1999), the ENSO index is defined as the 5-month running mean of deseasonalized 
SST anomalies averaged over the tropical eastern Pacific (defined as the area between 5°S-5°N 
and the South American coast and 180°W). The amplitude of interannual variation of LW 
broadband CRE is ~ 0.5 Wm" 2 . The running-mean time series of observed CRE and those of 
simulated by the AM2 and CanAM4 models are positively correlated with the ENSO index. They 
tend to have positive anomalies when the eastern Pacific SST anomalies are also positive (i.e., 
the El Nino state). GEOS-5 (red line in Fig. 8a) shows a different response to the ENSO index. 
The correlation between GEOS-5 broadband LW CRE time series and the ENSO index is slightly 
positive (~0.02) and but not statistically significant. The variation of each band CRE anomalies 
with the ENSO index closely follows the variation of LW broadband CRE anomalies, as shown in 
Fig. 8b for Bandl (results from other bands not shown here). 

When the 5-month running mean of the deseasonalized anomalies of f C RE is correlated 
with the ENSO index (Fig. 9), all three models agree with the observation on the signs of 
correlations over all bands excepts Band 4 (the ozone band). The correlation with the ENSO 
index is positive for Band 1 (the H 2 0 band) and negative for Bands 3 and 5 (the window bands). 
Correlations using the observations tend to be smaller than those from the simulations, which 
is likely due to the more noisy nature of the observations. The correlations in Fig. 9 indicate to a 
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large extent how the cloud top height varies at the interannual time scale with the ENSO index. 


Both the observation and models suggest that, averaged over the entire tropical ocean, the 
mean infrared-effective cloud top averaged over the tropical oceans tends to be elevated when 
the ENSO index is positive (El Nino state) and vice-versa. 

4. Model intercomparisons and comparisons with observations: spatial distributions 

Discussions in the previous section were about the averages over the entire tropical 
oceans. In this section, spatial distributions of band-by-band CRE and f C RE are to be discussed. 
As explained in Sectionl, the map of absolute CRE of a particular band largely resembles the 
map of broadband LW CRE (right panels in Fig. 2) and highlights the contrast of absolute CRE 
between regions featuring high clouds (where broadband LW CRE is ~ 60-80 Wm" 2 ) and regions 
dominated by low clouds (where broadband LW CRE is ~10-20 Wm" 2 or even less). Therefore, in 
this section the focus is on the spatial distributions of f C RE- Because the broadband and band-by- 
band CRE provide similar information, we mainly focus the discussion in this section to the 
spatial distributions of f C RE- 

Figure 10a shows the observed Bandl f C RE- As expected, the largest values of f C RE (~0.25- 
0.35) are found over the regions with frequent occurrence of high clouds, such as ITCZ, SPCZ, 
and Indian Monsoon regions. The smallest values occur in regions frequently covered by marine 
stratus (i.e. low cloud), such as the Pacific coast off South America, the Namibia coast, and the 
ocean region west of Australia. Figures lOb-lOd show the differences between modeled and 
observed Bandl f CRE . Simulated climatological mean positions of ITCZ over the equatorial Pacific 
and Atlantic are south to the observed ones, and the widths of ITCZ are also different, which 
accounts for the negative differences around 8°N and the slight positive differences between 0- 
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10 °N of the tropical Pacific and Atlantic. For regions where low clouds are prevalent, modeled 


f C RE are usually higher than the observed ones but the three models largely agree with each 
other. 

For Band2 (C0 2 band. Fig. lOe), the contrast in f CRE between high-cloud and low-cloud 
regions (0.16 vs. 0.11) is much smaller than that of Bandl (0.3 vs. 0.05). The models largely 
disagree with each other (Fig. 10 f-h) on the spatial distribution of f CRE . Compared to the 
observation, AM2 considerably underestimates Band2 f CRE over the low-cloud regions but 
largely agrees with observations over the rest of tropical oceans. GEOS-5 uniformly 
overestimates Band2 f CRE everywhere over the tropical oceans by ~0. 02-0. 04 with slightly larger 
overestimates for parts of the low-cloud regions. CanAM4 exceeds observations by ~0. 02-0. 04 
for the tropical oceans except the low-cloud region, where the difference is indeed slightly 
negative (~ -0.02 to 0). 

The observed Band3 f CRE (Fig. 11a) peaks over the low-cloud regions and drops over the 
high-cloud regions. This is consistent with the variation of Band3 f CRE with cloud top pressure 
(e.g.. Fig. 7a). Compared to the observation, the AM2 model overestimates Band3 f CRE nearly 
uniformly over the entire tropical oceans by ~0.04. The GEOS5 and CanAM4, on the other hand, 
underestimate Band3 f CRE over the low-cloud regions and are close to the observations for the 
rest regions of the domain with only slight positive or negative difference. These model- 
observation differences are opposite to those for Band2, for which the GFDL AM2 
underestimates over the low-cloud region and the GEOS5 and CanAM4 overestimate over 
majority of the tropical oceans. This is another way to expose the compensating biases among 
different bands in a GCM. As mentioned in Section 1 and Fluang et al. (2010), such differences 
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in modeled f C RE can be attributed to the differences in cloud top temperature. The differences 


between the AM2 and observations over both Band2 and Band3 consistently indicate that, on 
average, the cloud top temperature over the low-cloud regions are higher than observed ones, 
which leads to an overestimation of Band3 f CRE and the underestimation of Band2 f CRE . Similarly, 
for the GEOS-5, the mean cloud top temperature are lower than the observed ones, accounting 
for the underestimation of Band3 f CRE and the overestimation of Band2 f CRE . Note it is known 
that there is still difficulty in simulating temperature profiles in the low-cloud regions 
(temperature inversion, boundary layer mixing, etc.). So difference between simulated and 
observed cloud top temperature could be due to cloud top height difference as well as 
boundary-layer temperature difference between models and observations. For the CanAM4, 
the overestimation for Band2 is generally smaller than in GEOS-5 over the low-cloud regions 
but still noticeable over certain areas, such as the Peru coast and Namibia coasts, which 
suggests lower cloud top temperature than observed ones in these regions. 

The observed spatial map of Band4 f CRE (ozone band. Fig. lOe) is similar as that of Band3, 
owing to the fact that the ozone band is sensitive to the cloud and surface thermal contrast in a 
similar way as neighboring bands (i.e. Band3 and Band5). The ozone band in the AM2 model 
has a bandwidth 30 cm" 1 shorter than that used in the observations and the GEOS-5 and 
CanAM4 models (Table 2), which accounts for the consistently smaller contribution in the AM2 
(Fig. lOf) than in the observation. GEOS-5 agrees best with observations, even though smaller 
biases can still be seen over low-cloud regions consistent with those shown in Band3 over the 
same region. The calculated on-line ozone profiles continuously updated with the dynamical 
fields in the GEOS-5 simulation probably contribute to such good agreement. The spatial maps 
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of Band5 f C RE and model-observation difference (not shown here) are similar to those of Band3 


because the majority of the radiation of Band5 comes from the second window region (1100- 
1200 cm 1 ). 

As mentioned in Section 2, the GFDL AM2 radiation scheme assumes all clouds to be 
non-scattering in the LW while the NASA GEOS-5 and CCCma CanAM4 explicitly take the 
scattering into account (though the treatments are different in GEOS-5 and CanAM4). 
Therefore, it is meaningful to estimate how much of the difference shown in Figs. 10 and 11 is 
due to the inclusion of scattering. Similar to Fig. 7, we use the same conceptual model to 
estimate such difference. In one case cloud is assumed to be non-scattering and, in another 
case, the scattering is solved by a 4-stream DISORT solver in the MODTRAN5 package (Berk et 
al., 2005). f C RE of the non-scattering case is shown in Fig. 12a and the difference between the 
two cases (scattering - non-scattering) is shown in Fig. 12b. For clouds below 10km, the 
difference is negligible except for Band3 and cloud top below 4km. For clouds above 10km, 
there are noticeable differences in all five bands. From 10-13km, the difference is within ±0.03. 
The largest difference (~-0.09) is seen in Bandl for cloud top height at 14km. Such large 
negative difference can be understood from two aspects: (1) according to the Wien's 
displacement law, for a cloud top height at 14km (i.e., cloud top pressure at ~150mb and cloud 
top temperature around 210K), the fractional contribution from Bandl is as large as that from 
Band3 (Fig. 12a); (2) for ice, the imaginary part of the index of refraction of has a minimum 
value at 410cm 1 (Warren 1984; Warren & Brandt 2008), which means the inclusion of 
scattering would lead to the strongest scattering effect around this wavenumber. Overall, Fig. 
12b indicates that inclusion or not of scattering in the radiation scheme should have minor or 
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negligible impact on the f C RE, unless the cloud top is very high (>14km). Note over the tropical 


oceans, only a small fraction clouds can have their tops at or above 14km (~150 hPa). Cloud 
penetrating radar observation suggests that this fraction is 10% or less and is concentrated in 
the core of ITCZ area (Haynes & Stephens 2007). For the GCMs examined here, the monthly 
mean cloud fraction above 150 hPa is about the same as the observed. Therefore, the 
discrepancies between models and observations for either averages or spatial distributions, 
are primarily due to the difference in the modeled and observed cloud properties themselves, 
and not due to the inclusion of scattering or the details of scattering treatment in the radiation 
scheme. 

5. Conclusion and discussion 

Taking advantage of a multi-year multi-band LW CRE dataset over the tropical oceans 
derived from collocated AIRS and CERES observations, this study evaluates in details the LW 
band-by-band CRE simulated by three different GCMs. It quantitatively demonstrates that, even 
when the model-simulated broadband LW CRE largely agrees with the observed one, 
compensating biases from different LW bands can be prominent. The band-by-band CRE 
differences between models and observations can be as large as the difference in the LW 
broadband CRE. We prefer to use in our analysis the fractional contribution of each band to the 
LW CRE, f C RE, which has the unique advantage of being sensitive to cloud top 
height/temperature and not sensitive to cloud fraction. This property of f C RE enables us to 
separate the contribution of cloud fraction and cloud top height to the top-of-atmosphere LW 
flux. The averaged f CRE over the tropical oceans, their seasonal cycles, and their interannual 
variations are also examined. Bandl f C RE (H 2 0 band) is strongly anti-correlated with Band3 f C RE 
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and Band5 f C RE (window regions). The observed seasonal cycle is consistent with the seasonal 


cycles of ISCCP cloud top pressure climatology on both the phase and the amplitude. The 
simulated seasonal cycles of f CRE agree with observations on the phase but not on the 
amplitude which is much smaller in the simulatons, indicating smaller seasonal variations of 
simulated cloud top pressure, on average. As for interannual variations, observations and all 
GCMs agree on the positive (negative) correlations between Bandl f CRE (Band3 and Band5 f CRE ) 
and the ENSO index. Spatial distributions of band-by-band f CRE highlight the large discrepancies 
between models and observations in the regions with frequent occurrence of low clouds. While 
the total LW CRE is usually not used for the study of low clouds, fractional CRE contributions are 
sensitive to the low cloud top height (or equivalently, low cloud temperature). This 
complements the diagnostics of low cloud with shortwave reflectance, which is more sensitive 
to cloud liquid path than to the cloud top height/temperature. Cloud top height only directly 
affects the shortwave reflectance by changing the amount of solar absorption due to water 
vapor above the cloud. But change of cloud top height dominates the change of LW band-by- 
band f CRE . Therefore, such band-by-band f CRE offers a useful diagnostics for the low clouds 
simulated by the models. 

While this analysis focuses on fc RE for the reason articulated in Section 1, both absolute 
LW CRE and f CRE of each band are important to diagnose the quality of cloud simulations. Cloud 
fraction significantly affects the absolute LW CRE while f CRE is more sensitive to cloud top 
height/temperature. If a model is poor in producing correct cloud fractions for a certain type of 
cloud (e.g., the underestimation of low cloud fractions in the GEOS-5 model as noted in 
comparison with observed low cloud fractions, ), it would be logical and practical to first 
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address the issue of cloud fraction and then evaluate the f C RE- For both model evaluation and 


cloud feedback study, f C RE is a meaningful metric to use and it complements the absolute 
broadband LW CRE. 

To our knowledge, this is the first study to compare the simulated band-by-band LW CRE 
from multiple GCMs. The good agreement between NASA GEOS-5 and observation in the 
ozone-band CRE (Fig. 11) suggests that, when the detailed band-by-band decomposition of LW 
CRE are scrutinized against observations, the effect of ozone cannot be neglected. This has 
implication for the tuning of GCMs to reach the TOA energy balance as well as the climate 
projections in the presence of ozone recovery. Generally speaking, it is possible, if not very 
likely, to tune GCM to achieve a the TOA balance between outgoing longwave flux and net 
incoming shortwave flux for current climate while the band-by-band decomposition of LW flux 
(or LW CRE) is different from that observed. For example, the bias in the TOA energy imbalance 
due to simple representation of time-varying ozone fields can be "tuned" away by adjusting 
some less-constrained parameters in various GCM parameterization schemes. What is unknown, 
however, is what the consequence would be for future climate projection simulated by the 
same GCM when such incorrect band-by-band composition of TOA flux for the current climate 
exists. Our future work will try to address this kind of question by analyzing the band-by-band 
flux and CRE from simulations carried out for different future scenarios prescribed by the IPCC 
assessment. 
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Tables 


Table 1. Difference between OLR estimated from AIRS spectra with algorithms described in 
Huang et al. (2008) and Huang et al. (2010), denoted as OLR A | RS _huang, and OLR from collocated 
CERES measurements, denoted as OLRceres- Only observations over the tropical oceans are used 
here. Number of collocated observations as well as the mean difference and standard deviation 
are listed for each year as meanistandard deviation. 



Clear sky over the ocean 

Cloudy sky over the ocean 


Number of 
observations 

OLR A |RS_huang “ OLRceres 
(Wrrf 2 ) 

Number of 
observations 

OLR A iRs_huang ~ OLRceres 
(Wm' ! ) 

2003 

1.05x10 s 

0.84±1.56 

1.18xl0 7 

2.6615.68 

2004 

1.08x10 s 

0.6711.52 

1.24xl0 7 

2.1515.51 

2005 

1.08x10 s 

1.4211.68 

1.26xl0 7 

2.8815.79 

2006 

1.08x10 s 

1.5511.74 

1.27xl0 7 

2.9415.92 

2007 

1.04x10 s 

1.7311.86 

1.27xl0 7 

3.0015.92 


Table 2. The LW bands used in each GCM as well as the five bands used in this study for the 
comparisons across all GCMs and observations. The bandwidths are given in wavenumber (cm 
1 ). The major absorbers of each band are also listed. 


Major 

absorber(s) 

Band used in this 
study (cm' 1 ) 

GFDL Band ID 

G EOS-5 Band ID 

CanAM4 Band ID 

H 2 0 

Bandl: 0-560 &> 
1400 

Bands 1 & 8 
(0-560 & > 1400) 

Bandsl-2 (0-540) & 
Bands 8-9 (>1380) 

Bands 1-3 (>1400) 
Bands 8-9 (0-540) 

C0 2 

Band2: 560-800 

Band 2 
(560-800) 

Bands 3&10 
(540-800) 

Band 7 (540-800) 

H20 

continuum 

Band3: 800-980 

Band 3-4 
(800-990) 

Band4 

(800-980) 

Band 6 (800-980) 

o 3 

Band4: 980-1100 

Band 5 
(990-1070) 

Band5 

(980-1100) 

Band 5 (980-1100) 

H20 

continuum; 

ch 4 , n 2 o 

Band5: 1100-1400 

Band 6 (1070-1200) 
Band 7 (1200-1400) 

Band6 (1100-1215) 
Band7 (1215-1380) 

Band 4 (1100-1400) 
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Table 3. Multi-year means of LW broadband and band-by-band CRE averaged over the tropical 


oceans. The numbers in parentheses are the fractional contribution to the total LW broadband 
CRE (f CRE ). The AM2 results for Band4 are shaded and in italics to denote the prominent band- 
width differences between the GFDL AM2 model and the other two GCMs and observations 
(refer to Section 2 for detailed discussions). 



Observation 

GFDLAM2 

NASA GEOS-5 

CCCma CanAM4 


(Wm" 2 ) 

(Wm" 2 ) 

(Wm" 2 ) 

(Wm" 2 ) 

LW broadband 

28.5 (100%) 

28.0(100%) 

28.3 (100%) 

27.3 (100%) 

Bandl: 0-560 &> 1400 
cm 1 

5.59 (19.6%) 

5.31 (19.0%) 

5.11(18.1%) 

4.47 (16.3%) 

Band2: 560-800 cm" 1 

4.33 (15.2%) 

3.72 (13.3%) 

5.07 (18.0%) 

4.82 (17.6%) 

Band3: 800-980 cm" 1 

9.25 (32.4%) 

9.96 (35.6%) 

9.03 (32.0%) 

8.80 (32.2%) 

Band4: 980-1100 cm' 1 

3.72 (13.0%) 

1.68(6.0%) 

3.63 (12.9%) 

3.75 (13.7%) 

Band5: 1100-1400 cm" 1 

5.62 (19.7%) 

7.32 (26.2%) 

5.41 (19.2%) 

5.50 (20.1%) 
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Figure Captions 


Figure 1. (a) The percentage flux difference due to differences between intrinsic GCM 
bandwidths defined in each GCM and the bandwidths used in this study (refer to Table 2 for the 
details). Flux is calculated for a blackbody at 220K. (b) Same as (a), but for a blackbody at 298K. 
(c) The flux difference between (a) and (b) expressed in percentage. 

Figure 2. The multi-year mean of all-sky OLR (left panels) and LW broadband LW CRE (right 
panels) derived from collocated AIRS and CERES observations, the GFDL AM2, NASA GE0S5, and 
CCCma CanAM4 model simulations. Note that color scales are different for the left panel and 
the right panels. 

Figure 3. The mean seasonal cycles of Bandl CRE (blue lines), Band2 CRE (green lines), and 
Band3 CRE(red lines). Definitions of bandwidths can be found in Table 2. For better visualization, 
CREs from the observation and GEOS-5 model are plotted on one scale (±lWm~ 2 ) and those 
from the AM2 and CanAM4 model are plotted on another scale (±0.5Wm 2 ). 

Figure 4. The mean seasonal cycle of f C RE, the fractional contribution of a band to the total LW 
CRE. Band 1-3 are plotted in blue lines, green lines, and red lines, respectively. 

Figure 5. Correlation coefficients between the mean seasonal cycle of Bandl and that of all 
other bands. The observed relation is plotted in blue dotted line. The GFDM AM2, NASA GEOS-5, 
and CanAM4 models are plotted in green, red, and cyan dotted lines, respectively. 

Figure 6. (a) The ISCCP climatological seasonal cycle (deviations from the climatological mean) 
of cloud fraction for high cloud (blue), middle cloud (green), and low cloud (red) averaged over 
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the tropical oceans, (b) The ISCCP climatological seasonal cycle of cloud top pressure averaged 


over the tropical oceans. 

Figure 7. (a) Change of Bandl f CRE and Band3 f CRE with respect to the cloud top pressure in a 
simple conceptual model in which cloud is assumed to be optically thick and typical tropical 
sounding profiles are used. Radiative calculation is done with MODTRAN5. (b) Observed 
standard deviations of the seasonal cycles of Bandl f CRE and Band3 f CRE vs. estimated 
counterparts. Black dash line denotes the 1:1 slope. 

Figure 8. (a) 5-month running means of deseasonalized anomalies of total LW CRE from the 
observations (blue line), the GFDL AM2 (green line), NASA GEOS-5 output (red line), and CCCma 
CanAM4 (cyan line) simulations. The ENSO index is plotted in gray dash line, (b) Same as (a) 
except for the Bandl CRE instead of total LW CRE. Following previous studies, the ENSO index is 
defined as 5-month running mean of the SST anomalies averaged over the region of 5°S-5°N, 
the South American coast and 180°W. 

Figure 9. Correlation coefficients between the ENSO index and the 5-month running means of 
deseasonalized f CRE of each band. Results from observation and three models are shown in 
different colors. Figure 10 (a) Map of observed long-term mean of Bandl f CRE . (b)-(d) Maps of 
model-observation difference in the long-term mean of Bandl f CRE for the GFDL AM2, NASA 
GEOS-5, and CCCma CanAM4 simulations, respectively. Note the same color scale is used for 
(b)-(d). (e) same as (a) except for Band2 f CRE . (f)-(h) same as (b)-(d) except for Band2 f CRE . 
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Figure 11. Same as Figure 10 except for Band3 f C RE and Band4 f C RE- 


Figure 12. (a) Band-by-band f C RE for different cloud top height (2km-14km). Clouds are assumed 
to be opaque and non-scattering, (b) The difference in f CRE between the case of scattering 
cloud and the case of non-scattering cloud when everything else remains the same. 
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T = 220K 



Figure 1. (a) The percentage flux difference due to differences between intrinsic GCM 
bandwidths defined in each GCM and the bandwidths used in this study (refer to Table 2 for the 
details). Flux is calculated for a blackbody at 220K. (b) Same as (a), but for a blackbody at 298K. 
(c) The flux difference between (a) and (b) expressed in percentage. 
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Figure 2. The multi-year mean of all-sky OLR (left panels) and LW broadband LW CRE (right 
panels) derived from collocated AIRS and CERES observations, the GFDL AM2, NASA GEOS5, and 
CCCma CanAM4 model simulations. Note that color scales are different for the left panel and 
the right panels. 
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Figure 3. The mean seasonal cycles of Bandl CRE (blue lines), Band2 CRE (green lines), and 
Band3 CRE(red lines). Definitions of bandwidths can be found in Table 2. For better visualization, 
CREs from the observation and GEOS-5 model are plotted on one scale (±lWm 2 ) and those 
from the AM2 and CanAM4 model are plotted on another scale (±0.5Wm 2 ). 
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Figure 4. The mean seasonal cycle of f C RE, the fractional contribution of a band to the total LW 
CRE. Band 1-3 are plotted in blue lines, green lines, and red lines, respectively. 
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Figure 5. Correlation coefficients between the mean seasonal cycle of Bandl and that of all 
other bands. The observed relation is plotted in blue dotted line. The GFDM AM2, NASA GEOS-5, 
and CanAM4 models are plotted in green, red, and cyan dotted lines, respectively. 
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Figure 6. (a) The ISCCP climatological seasonal cycle (deviations from the climatological mean) 
of cloud fraction for high cloud (blue), middle cloud (green), and low cloud (red) averaged over 
the tropical oceans, (b) The ISCCP climatological seasonal cycle of cloud top pressure averaged 
over the tropical oceans. 


38 




Figure 7. (a) Change of Bandl f C RE and Band3 f C RE with respect to the cloud top pressure in a 
simple conceptual model in which cloud is assumed to be optically thick and typical tropical 
sounding profiles are used. Radiative calculation is done with MODTRAN5. (b) Observed 
standard deviations of the seasonal cycles of Bandl f C RE and Band3 f C RE vs. estimated 
counterparts. Black dash line denotes the 1:1 slope. 
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Figure 8. (a) 5-month running means of deseasonalized anomalies of total LW CRE from the 
observation (blue line), the GFDL AM2 (green line), NASA GEOS-5 output (red line), and CCCma 
CanAM4 (cyan line) simulations. The ENSO index is plotted in gray dash line, (b) Same as (a) 
except for the Bandl CRE instead of total LW CRE. Following previous studies, the ENSO index is 
defined as 5-month running mean of the SST anomalies averaged over the region of 5°S-5°N, 
the South American coast and 180°W. 
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Figure 9. Correlation coefficients between the ENSO index and the 5-month running means of 
deseasonalized f C R E of each band. Results from observation and three models are shown in 
different colors. 
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Figure 10 (a) Map of observed long-term mean of Bandl f C RE- (b)-(d) Maps of model-observation 
difference in the long-term mean of Bandl f C RE for the GFDL AM2, NASA GEOS-5, and CCCma 
CanAM4 simulations, respectively. Note the same color scale is used for (b)-(d). (e) same as (a) 
except for Band2 f C RE- (f)-(h) same as (b)-(d) except for Band2 f C RE- 
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Figure 11. Same as Figure 10 except for Band3 f C RE and Band4 f C RE- 
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Figure 12. (a) Band-by-band f CRE for different cloud top height (2km-14km). Clouds are assumed 
to be opaque and non-scattering, (b) The difference in f CR E between the case of scattering 
cloud and the case of non-scattering cloud when everything else remains the same. 
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